Quasi-one-dimensional magnons in an intermetallic marcasite 
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We present inelastic neutron scattering measurements and first principles calculations examining 
the intermetallic marcasite CrSb2. The observed spin wave dispersion implies that the magnetic 
interactions are strongly one-dimensional with antiferromagnetic chains parallel to the crystalline 
c-axis. Such low-dimensional excitations are unexpected in a semiconducting intermetallic system. 
Moreover this material may be further interesting in that the magnetic anisotropy may enhance 
thermoelectric properties along particular crystallographic directions. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.30.Et, 75.30.Ds, 78.70.Nx 



Efficient and cost-effective power generation has re- 
vived interest in thermoelectric materials in an energy 
resource dominated economy. CrSb2 has gained re- 
cent attention as a narrow gap semiconductor with po- 
tential thermoelectric properties. Investigations into 
substitutional compounds such as CrSb2_xTe2; P, [^j, 
CrSb2-a;Sna;[3l and Cri_a;Rua;Sb2 [1] have found changes 
in the electronic properties, with an increase in the ther- 
moelectric figure of merit for Te substitution. CrSb2 has 
also been considered as an electrode material for lithium 
based batteries showing improvement compared to us- 
ing pure antimony. ^ Q Here we report inelastic neutron 
scattering (INS) measurements and first principles cal- 
culations for this compound. We find that CrSb2 has 
significant exchange anisotropy exhibiting a quasi-one- 
dimensional (Id) magnon spectrum. This spectrum has 
implications for the thermal conductivity as a function of 
crystallographic direction. In addition, the determined 
exchange constants and anisotropy allow one to classify 
CrSb2 as an S' = 1 antiferromagnetic (AFM) chain with 
interchain coupling and on-site anisotropy. 

CrSb2 is an intermetallic marcasite where the Cr atoms 
are located at the body center and corner positions of 
the orthorhombic unit cell with room temperature lat- 
tice constants a = 6.028, b = 6.874, and c = 3.272 1.0 
The Cr atoms are at the center of a distorted octa- 
hedron of Sb sites as shown in Fig. [T] Early ther- 
modynamic measurements found CrSb2 to have AFM 
correlations!^ Q, but the ordered magnetic phase with 
T/v = 273(2) K was not observed until much later. Q The 
ordered magnetic structure is shown in Fig. [T] where the 
moments are roughly perpendicular to the (101) plane. 
This magnetic structure results in magnetic Bragg peaks 
at the {H ^^) positions where H, K, and L 

are the Miller indices of the chemical unit cell. The 
body centered Cr sites form the identical AFM struc- 
ture as the corner sites, albeit displaced along the [111] 
direction. The magnetic unit cell is doubled along the 
b and c axes with an ordered moment corresponding to 



spin 5* = 0.97(2) as determined from neutron diffrac- 
tion measurements. Q Early band structure considera- 
tions have shown that there are two localized d electrons 
per Cr'*+ in CrSbz. [l^ 

Local spin density approximation (LSDA) calculations 
were performed using the general potential linearized 
augmented planewave (LAPW) method. [ll| Local orbital 
extensions were employed to include the high lying semi- 
core states and to relax any linearization errors. The cal- 
culations included no shape approximations to either the 
potential or charge density. Relativity was included for 
the valence states within a scalar relativistic approxima- 
tion, while full relativity was included for the core states, 
within an atomic like approximation. Well converged ba- 
sis sets consisting of more than 850 LAPW functions 
plus local orbitals were used for the primitive unit cell 
(two formula units) of CrSb2. Brillouin zone samplings 
were done with the special k-points method, using 128 




FIG. 1: (color online) Nuclear and magnetic structure of 
CrSb2.[3, 01 Red and green spheres are Cr sites, where the 
vector and color represent the spin orientation. Small blue 
spheres are the Sb sites shown only for the three unit cells on 
the right of the figure for clarity. Nearest neighbor exchange 
constants are labeled and shown as heavy black lines. 
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k-points in the irreducible wedge of the primitive cell and 
equivalent samplings for the larger AFM cells. 

We consider ferromagnetic (FM) and several AFM 
ground states including the ordered magnetic structure 
found for CrSb2, which we denote as AFM-D. Other 
AFM structures include doubled magnetic cells along the 
a (AFM-A) and h (AFM-B) axes, and FM correlations 
within each of the sublattices but AFM relative to one 
another (AFM-G). The lowest calculated LSDA energy 
is found for AFM-D, -0.522 eV per formula unit relative 
to the non-spin polarized case, with other values listed 
in Fig. [51 All of the configurations were found to be 
conducting except for AFM-D. The calculated density of 
states for AFM-D is shown in Fig. The calculated 
bandgap of 0.07 eV agrees with low temperature resis- 
tivity measurements. [l2| The large energy difference be- 
tween the AFM-D and FM structures compared to the 
relatively low T^r in a lattice with no geometrical frus- 
tration suggests interesting magnetic correlations. The 
energy of the AFM-G state is 0.025 eV per Cr higher 
than the FM state, indicating that the inter-sublattice 
coupling, Jill in Fig.[Tl is very weak and likely FM. The 
energy of the other AFM states implies weak FM ex- 
change along a, almost no interaction along 6, and strong 
AFM exchange along c in agreement with the ordered 
magnetic structure. The large AFM exchange along the 
c axis would make the magnetic correlations in CrSb? 
quasi-ld, thus explaining the suppressed value of Tjv.[l3| 
This also explains the difficulty in observing the phase 
transition in prior thermodynamic measurements. The 
Id magnetism would result in much of the magnetic en- 
tropy being exhausted with a very broad hump in sus- 
ceptibility as the local order develops, with only smaller 
features at T^f. 

Single crystal INS provides an efBcient probe of the 
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FIG. 2: LSDA calculated electronic density of states N{E) per 
formula unit (fu) for the AFM order of the magnetic ground 
state (AFM-D). The table inset lists the configurations and 
calculated LSDA energies per formula unit relative to the non- 
spin-polarized case (see text). 



LSDA predictions of the exchange coupling for CrSb2. 
Single crystals were grown from an Sb flux with a start- 
ing composition of CrSbig. The material was heated in 
vacuum in a 10 cc alumina crucible to 1273 K, held for 
24 h, and then cooled to 913 K at 1 K/hr. The excess Sb 
liquid was centrifuged away, leaving crystals with typical 
dimensions of 5x5x5 mm'^. Single crystal INS measure- 
ments of CrSb2 were performed at the HB3 triple-axis 
spectrometer (TAS) and the ARCS chopper spectrome- 
ter at Oak Ridge National Laboratory's High Flux Iso- 
tope Reactor and Spallation Neutron Source respectively. 
Measurements at IIB3 were performed for samples in the 
(0 K L)[m « 4.05 g] and (H L L)[m « 3.5 g] scatter- 
ing planes. The instrument was configured with a verti- 
cally focused Be (002) monocrhomator and a PG (002) 
analyzer with fixed final energy Ef = 30.5. Pyrolytic 
graphite filters were in place between the sample and ana- 
lyzer. Constant wave vector, Q, and energy transfer, Uuj, 
scans were performed primarily along the (H k/2 k/2), 
(0 K k/2) and (0 k/2 L) directions, where k = 1 or 3 be- 
tween 5 and 70 mcV. Data were acquired for a fixed mon- 
itor count unit (mcu) to account for variations in neutron 
flux as a function of incident energy, Ei . Measurements 
at ARCS were performed with the (0 K L) plane in the 
equatorial plane of the instrument with Ei = 150 meV 
neutrons. Data were measured by rotating about the 
vertical (10 0) axis in increments of 0.75 degrees over a 
total of 62 degrees with the incident wave vector initially 
pointed 9 degrees from the (0 1 0) direction. Data were 
normalized to the accumulated charge on the source and 
histogramed into Q and hto space. 

Figure [3] shows the measured scattering intensity as a 
function of Tilu and Q for CrSb2 at T = 10 K. Given the 
inversion symmetry of the magnetic structure, the data 
from ARCS, Fig. [3l^d)-(f), have been averaged for data 
acquired at ±Q and for higher Q Brillouin zones scaled 
by the Cr magnetic form factor. Data have been back- 
ground subtracted using data scaled from larger Brillouin 
zones to suppress phonon scattering. There is a clear 
single excitation that disperses along all three principal 
axes. The overall minimum in the dispersion appears at 
w 25 meV at (i/, ^i^, ^L±l) for integer H,K and L. 
The kinematic limits of the TAS measurements did not 
allow us to follow this mode above 60 meV for small val- 
ues of |Q|. The ARCS measurements are able to observe 
this mode up to its maximum energy of w 78 meV. The 
bandwidth along the H and K direction is 13.3 and 14.8 
meV respectively, and the bandwidth of the excitation 
along L is 53.6 meV. 

Figure m shows the (OKL) scattering intensity from the 
ARCS measurement integrated for 25 < hw < 100 meV. 
The scattering intensity falls off rapidly with | Q | indicat- 
ing magnetic fluctuations. There are ridges of scatter- 
ing consistent with the greatest dispersion being along 
the c-axis as shown in Fig. |31 Prior investigations of 
the high temperature magnetic susceptibility employed 
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FIG. 3: (Color online) T = 10 K INS spectra of CrSb2 along principal axes of the reciprocal lattice, (a)-(c) [(d)-(f)] Data were 
measured using IIB3[ARCS]. (d)-(f) Have been background subtracted as described in the text. (a)[(b)-(f)] were measured with 
the sample in the (H K K)[(0 K L)] scattering plane. Open circles correspond to peak positions determined from extracted 
constant Q scans as described in the text. Open squares in (c) correspond to peak positions from ARCS and HB3 constant Q 
scans along the {H ^^"'^^'"^ L) directions where H and K are integers. VerticalfHorizontal] error bars represent the fitted peaks 
full width at half maximum [range of Q integration]. The grey solid line is the fitted spin- wave dispersion described in the text. 



a two-state magnetic system with an energy gap of 70.7 
meV to describe the results. [13] Although this value is 
within the energy scale of the measured magnon, it is 
quite far from the peak in the magnetic density of states. 
Our INS measurements provide a picture of the micro- 
scopic interactions in this material, and clearly shows 
that the gap in the magnetic spectrum is approximately 
25 meV. 

We project the data from ARCS and HB3 into a col- 
lection of constant Q scans by integrating over ±0.05 rlu 
in each of the principal axes directions. These individual 
scans are then fit to a single inelastic Gaussian peak for 
20 < ?j,a; < 80 meV, allowing us to extract the dispersion 
of the magnetic excitation. A portion of these fitted peak 
positions is shown in Fig. [3l The points in Fig. [Sjc) arc 
from a series of Q's along the {r] '^'^2^ L) direction where 
77 and K are integer values measured with both HB3 and 
ARCS. This panel illustrates the extent of the dispersion 
along the L direction. Although the electronic band-gap 
is small, we consider the cooperative magnetism in CrSb2 
to be local moment in nature. This allows one to use lin- 
ear spin wave theory (LSWT) to calculate a dispersion 
relation, huj(Q)sw for comparison to the measured data. 



The tilts and distortion of the Cr-Sb octahedron, the di- 
rection of the ordered moments, and the gap in the AFM 
spectrum implies the need for an anisotropy parameter, 
D. The Heisenberg Hamiltonian then becomes 



(1) 



where a is summed over a, h, and c, and we define the 



axis to be along the ordered moment. [1J,[15| This results 
in the spin-wave dispersion 



= 2S[Ja[cos{2TiH)-\\ + .h + Jc + D] 
Bq = 2S[JbCOs{2'!TK) + JcCos{2ttL)] 



(2) 



where S is the magnitude of the spin quantum and the 
exchange constants are illustrated in Fig. [1] Fitting to 
huj{Q)sw with S = 1 results in the solid grey line in 
Fig. [3] and the parameters Jc = 35.7(2), Ja = —1.31(6), 
Jb = 1.69(7) and D = 2.07(9) meV. Including a Jm 
exchange constant has negligible improvement on the 
of the fit. LSWT accounts well for the measured low- 
temperature dispersion in CrSb2. The resulting values 
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FIG. 4: (color online)Intensity in the (OKL) plane integrated 
for 25 < Tiuj < 100 meV and \H\ < 0.1 rlu (reciprocal lattice 
units). Data are from the ARCS measurements at T = 10 K. 
Solid lines are the contours of the calculated spin-wave scat- 
tering intensity as described in the text. Contour lines are 
plotted in logarithmic steps. White area in bottom left was 
not sampled during the measurement. 



also agree well with the FM and AFM spin-spin interac- 
tions based upon the low temperature ordered magnetie 
structure, and the predictions based upon LSDA calcu- 
lations discussed earlier. The Jc exchange is an order of 
magnitude larger than along the other axes in support of 
the predicted quasi- Id magnon behavior. 

We also examine the integrated scattering intensity in 
the {OKL) plane as shown in Fig. |4l We calculate the en- 
ergy resolution convolved one-magnon cross-section 14 1 
based upon the determined exchange constants and the 
resolution function of the ARCS spectrometer. We in- 
clude the Cr magnetic form factor, and plot the cross 
section as a contour in Fig. |4l The calculated scattering 
intensity agrees very well with the measured spectrum. 

For T > Tn, the magnetic excitations in CrSb2 could 
be considered 5 = 1 quasi-ld magnons propagating 
out of a spin-liquid ground state suggesting comparison 
to a Haldane chain dispersion. Thermal population of 
phonons makes an accurate comparison using INS diffi- 
cult at high temperatures, but the extracted exchange 
constants based upon LSWT for T < T/y place CrSb 



well within an ordered phase of the phase diagram. [16lll7 | 
Measurements of the narrow-gap semiconducting marc- 
asite FeSb2 have found this system to be nearly mag- 
netic, and that substitution with Cr as Fei_a;Cr2;Sb2 with 
X > 0.45 pushes the system into the AFM phase with 



Tn increasing with a;.[ii,[ig| Although INS of FeSb2 ob- 
served no clear magnetic excitations, [i^ investigations 
into suppressing the T/y of CrSb2 via chemical substitu- 
tion are fundamentally interesting to try to move through 
the Haldane phase diagram and ultimately access the 
spin-liquid Haldane phase. Another route to adjusting 
the exchange constants along with the electrical con- 



ductivity would be to change the atomic percent of Sb. 
The high temperature trend of the magnetic susceptibil- 
ity is very similar for the range of compositions between 
CrSbi.75 to CrSbs.i 

The quasi-ld magnetic excitation spectrum of CrSb2 
will affect its thermoelectric properties. Thermal trans- 
port measurements have been performed for other low- 
d magnets. In the case of the S = ^ spin-ladder 
Sri4_.TCaxCu?404l[2l| and the 5 = 1 Haldane chain 
NENP[22|, it was found that thermal transport was en- 
hanced along the Id axis by a factor of up to w 15. 
The large magnetic gap and spin-wave velocity found 
for CrSb2 would serve to enhance the thermal conduc- 
tivity along the Id axis for temperatures of the order of 
Jc- This implies that the room temperature thermoelec- 
tric figure of merit may be improved by orienting the Id 
axis orthogonal to the flow of charge. While the lattice 
thermal conductivity is probably dominant, the magnetic 
contribution is potentially important and it along with 
its anisotropy should be considered when assessing the 
thermoelectric properties of CrSb2. 
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